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ABSTRACT 

The non-axisymmetric features observed in the discs of dwarf novae in outburst are 
usually considered to be spiral shocks, which are the non-linear relatives of tidally 
excited waves. This interpretation suffers from a number of problems. For example, the 
natural site of wave excitation lies outside the Roche lobe, the disc must be especially 
hot, and most treatments of wave propagation do not take into account the vertical 
structure of the disc. 

In this paper I construct a detailed semi-analytical model of the non-linear tidal 
distortion of a thin, three-dimensional accretion disc by a binary companion on a 
circular orbit. The analysis presented here allows for vertical motion and radiative 
energy transport, and introduces a simple model for the turbulent magnetic stress. 
The Til — 2 inner vertical resonance has an important influence on the amplitude 
and phase of the tidal distortion. I show that the observed patterns find a natural 
explanation if the emission is associated with the tidally thickened sectors of the outer 
disc, which may be irradiated from the centre. According to this hypothesis, it may be 
possible to constrain the physical parameters of the disc through future observations. 
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^ ; 1 INTRODUCTION 



1.1 Tidal effects on accretion discs 

Accretion discs are commonly found in interacting binary 
star systems where matter is transferred from a normal star 
towards a compact companion (e.g. Lubow & Shu 1975). 
Tidal forces have an important role in these systems. The 
binary orbit is usually made accurately circular by the tidal 
interaction between the stars. The disc surrounding the com- 
pact star also experiences a strong tidal force from the com- 
panion, resulting in a non-axisymmetric distortion of the 
disc. Angular momentum transported outwards through the 
disc by 'viscous' stresses is transferred to the binary orbit 
through tidal torquesj^] 

The tidal interaction between a disc and an orbiting 
companion is also important in young binary stars and pro- 
toplanetary systems (e.g. Lin & Papaloizou 1993). In these 
situations there may be discs both interior and exterior to 
the companion's orbit, which need not be circular. 

Early theoretical studies addressed the tidal influence 
of a companion with a circular orbit and a mass ratio 
q — M% /Mi of order unity on a circumstellar disc (Paczyhski 



1 In reality the 'viscous' stress is most likely to be a turbulent 
magnetic stress resulting from the non-linear development of the 
magnetorotational instability (Balbus & Hawley 1998). 



1977; Papaloizou & Pringle 1977) . The method of Paczyhski 
(1977) makes use of the idea that, in a thin accretion disc, 
the effects of pressure and viscosity are relatively weak. In 
the absence of resonances, and provided that the trajecto- 
ries do not intersect, the motion of the gas is essentially 
ballistic. The ballistic trajectories in a circular binary sys- 
tem are just the orbits of the restricted three-body problem, 
which has been studied in great detail in celestial mechanics 
(e.g. Szebehely 1967). Paczyhski (1977) showed that there 
is a family of these orbits that lie in the binary plane and 
are periodic in the binary frame, enclosing the compact star 
(Fig. 1). They reduce to circular Keplerian orbits when the 
companion star is removed. It is therefore natural to suppose 
that the streamlines of a tidally distorted disc correspond 
accurately to Paczyhski's orbits. The maximum size of the 
disc can be identified as the largest non-intersecting orbit. 

An alternative way of analysing the tidal distortion of 
a disc is through linear perturbation theory. A solution is 
first obtained for an axisymmetric disc around an isolated 
star. The tidal force is introduced as a perturbation and 
the linearized equations are solved to determine the tidal 
distortion. A solution can be found that is steady in the bi- 
nary frame. The distortion of the disc is very slightly out 
of phase with the tidal forcing owing to the effective viscos- 
ity of the disc. This results in a tidal torque that increases 
with increasing radius until it becomes comparable to the 
local viscous torque and the disc is truncated. Such a cal- 
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Figure 1. Paczynski's orbits in the restricted three-body problem for a mass ratio 5 = 1. The prograde circumstellar case is depicted. 



culation was carried out for a two-dimensional disc model 
by Papaloizou & Pringle (1977). At least for mass ratios of 
order unity, the tidal truncation radii determined by the two 
methods agree reasonably well. 

1.2 Resonances 

Additional complications arise if the disc contains radii at 
which the tidal forcing resonates with a natural mode of 
oscillation of the disc. For mass ratios g<0.25, a band of 
Paczynski's orbits around the 3:1 resonance are parametri- 
cally destabilized (Hirose & Osaki 1990). In a continuous 
disc this causes a local growth of eccentricity (Lubow 1991), 
which may be able to compete with dissipative processes to 
sustain a global eccentric mode of the disc, as observed in 
superhump binaries. 

When the mass of the companion is much smaller, one 
or more Lindblad resonances may be reached before the disc 
is truncated. Such resonances usually result in the launching 
of non-axisymmetric waves (Goldreich & Tremaine 1979). 
Unless the waves are able to propagate a significant distance 
through the disc, the torque exerted at a Lindblad resonance 
is localized and independent of the viscosity. Such a localized 
torque will often truncate the disc in the case of a very low- 
mass companion (e.g. Lin & Papaloizou 1986). 

1.3 Three-dimensional effects 

The methods of Paczynski (1977) and Papaloizou & Pringle 
(1977) both neglect the three-dimensional nature of the disc. 
In a linear analysis, Lubow (1981) has shown that the verti- 
cal component of the tidal force induces homogeneous verti- 
cal expansions and contractions of the disc, which are char- 
acterized by a vertical velocity proportional to the distance 
z from the mid-plane. The natural frequency of this mode 
of oscillation, which corresponds to the n = 1 p mode in 
the notation of Lubow & Pringle (1993), or the p\ mode in 



the classification of Ogilvie (1998), depends on the adiabatic 
exponent 7. 

As in the case of horizontal forcing, the response to ver- 
tical tidal forcing may involve resonant and non-resonant 
components. The equivalent of a Lindblad resonance for the 
Pi mode is a vertical resonance^], the location of which de- 
pends on 7 (Lubow 1981). If such a resonance lies within the 
disc, a resonant torque is exerted there and waves may be 
excited. Vertical resonances in a thin disc in a binary sys- 
tem are intrinsically much weaker than Lindblad resonances. 
However, for mass ratios of order unity, the m = 2 inner 
vertical resonance typically lies within the tidal truncation 
radius, while the Lindblad resonances all lie well outside the 
disc (outside the Roche lobe, in fact). The vertical resonance 
can therefore have an important influence on the disc. 

There is also a non-resonant, non-wavelike response to 
vertical tidal forcing at all radii. Dissipation can result in a 
significant phase shift in this response, causing an additional 
non-resonant tidal torque. 

Although these effects have been widely neglected in the 
literature, Stehle & Spruit (1999) have recently proposed 
an interesting method of including vertical motion, in an 
approximate way, in two-dimensional numerical simulations 
of accretion discs. Stehle & Spruit (1999) also drew attention 
to the resonant excitation of vertical oscillations by the m = 
2 component of the tidal potential. 



1.4 Non-axisymmetric structure in accretion discs 

In recent years, non-axisymmetric features have been de- 
tected in the discs of several dwarf novae in outburst us- 
ing Doppler tomography (Steeghs, Harlaftis & Home 1997; 



2 Distinct from the vertical resonances that excite bending modes 
(e.g. Shu, Cuzzi & Lissauer 1983). Here we are concerned with 
motions that are symmetric about the mid-plane. 
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Steeghs 2001 and references therein). The tomograms indi- 
cate the surface brightness of the disc, in the selected emis- 
sion line, in a two-dimensional velocity space. The pattern 
generally appears as two crescents forming a broken ring, 
and the same is true when the pattern is mapped into real 
space, if circular Keplerian motion is assumed. 

Steeghs et al. (1997) interpreted the non-axisymmetric 
features as evidence of 'spiral structure', 'spiral waves' or 
'spiral shocks' in the disc. Observers of similar features 
in other systems have confidently reported them as spiral 
shocks (Joergens, Spruit & Rutten 2000; Groot 2001). Nev- 
ertheless, owing to the substantial complications in predict- 
ing line emission from cataclysmic variable discs (e.g. Robin- 
son, Marsh & Smak 1993), these identifications should be 
regarded with some caution. The tomograms may not accu- 
rately reflect the dynamics of the disc, because they originate 
in the tenuous layers of the atmosphere. It might even be 
questioned whether the observed features indeed resemble 
either a spiral or a wave. 

Two-armed spiral shocks had previously been seen in 
two-dimensional numerical simulations of the formation of 
very hot accretion discs by Roche-lobe overflow (Sawada, 
Matsuda & Hachisu 1986). Spruit (1987) showed with an 
analytical model that they could allow accretion to occur 
in an essentially inviscid disc, although the accretion rate is 
very small in a thin disc. Spiral shocks in two-dimensional 
disc models can be understood as the non-linear relatives 
of linear acoustic-inertial waves (Larson 1990). A popular 
conception is that such waves are excited by tidal forcing 
at the outer edge of the disc and steepen into shocks as 
they propagate inwards. Although the m = 2 inner Lind- 
blad resonance lies outside the Roche lobe, it is sometimes 
considered to launch waves through its 'virtual' effect on the 
disc (Savonije, Papaloizou & Lin 1994). 

However, an analysis of Lindblad resonances in three- 
dimensional discs (Lubow & Ogilvie 1998) indicates that the 
dominant wave that is launched is the fundamental mode of 
even symmetry, or f e mode. Only in a vertically isother- 
mal disc does this mode behave like the waves in a two- 
dimensional disc model; in a thermally stratified disc the f e 
mode behaves like a surface gravity mode and has quite dif- 
ferent propagational properties. Typically, it dissipates after 
travelling only a short distance from the Lindblad resonance 
(Lubow & Ogilvie 1998; Bate et al. 2001). 

Moreover, as pointed out by Godon, Livio & Lubow 
(1998), spiral waves in a thin accretion disc tend to be tightly 
wound except near resonances. If the non-axisymmetric fea- 
tures are to be interpreted as a two-armed wave, then the 
observations indicate at most one radial wavelength. Either 
the wave is highly localized in the radial direction, or it is 
not tightly wound. As the Lindblad resonances lie well out- 
side the Roche lobe in IP Peg, it is very difficult to account 
for unwound waves in two-dimensional disc models unless 
the disc is unrealistically hot. 

1.5 Vertical resonances 

However, as mentioned by Godon et al. (1998), a vertical 
resonance could lie inside the disc in IP Peg. In fact, the lo- 
cation of the m = 2 inner vertical resonance for an adiabatic 
exponent 7 = 5/3 agrees well with the observed features. In 
an inviscid, linear theory a pf-mode wave is emitted from 



the resonance (Lubow 1981), which would not be tightly 
wound near the site of launching. However, in a separate 
paper (Ogilvie, in preparation) I make a detailed analysis 
of this resonance allowing for non-linearity and dissipation, 
and show that the resonance is typically saturated, broad- 
ened and damped to the extent that the emitted wave can 
be neglected. Therefore no such wave will be considered in 
the present paper. Nevertheless, the presence of the verti- 
cal resonance still has an important role in determining the 
tidal distortion, even if a wave is not generated. 



1.6 Plan of the paper 

The purpose of this paper is to determine the non-linear 
tidal distortion of a thin, three-dimensional accretion disc by 
a binary companion on a circular orbit. The fluid dynamical 
equations are solved in three dimensions by means of asymp- 
totic expansions, using only the fact that the disc is thin. The 
resulting model is more complete than any previous analy- 
sis. Regarding the two-dimensional aspects of the problem, 
the model of Papaloizou & Pringle (1977) is extended by 
making it non- linear, while the model of Paczyhski (1977) is 
developed by including the effects of pressure and turbulent 
stresses. In addition, allowance is made for vertical motion, 
radiative energy transport and stress relaxation. That is, 
a Maxwellian viscoelastic model of the turbulent stress is 
adopted, which accounts for the non-zero relaxation time of 
the turbulence. 

This approach is part of the programme of 'continuum 
celestial mechanics' that has already addressed the dynam- 
ics of warped discs (Ogilvie 1999, 2000) and eccentric discs 
(Ogilvie 2001) using a closely related method. It provides 
a completely independent alternative to direct numerical 
simulations. The present method is more naturally suited 
to thin, but fully vertically resolved, discs. It also allows a 
more rapid solution, and is in some respects more advanced 
in physical content, than present-day simulations. 

Although the finer details of the solutions obtained are 
not at present observationally testable, comparisons can be 
made at a simple level between the non-axisymmetric struc- 
ture in the computed models and the features observed in 
Doppler tomograms. It will be seen that the observed pat- 
terns find a natural explanation if the emission is associated 
with the tidally thickened sectors of the outer disc, which 
may be irradiated from the centre. 

The remainder of this paper is organized as follows. In 
Section 2 Paczyhski's orbits are used to define a orbital co- 
ordinate system that is used in solving the fluid dynamical 
equations. In Section 3 the model for the turbulent stress is 
explained. The basic equations are written in the orbital co- 
ordinate system in Section 4. The asymptotic development 
and solution of the equations is worked out in Section 5, 
leading to a prediction for the non-axisymmetric distortion 
of the disc. Simple synthetic Doppler tomograms are com- 
puted and compared with observations in Section 6. A sum- 
mary and discussion of the results are found in Section 7. 
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2 ORBITAL COORDINATES 

Consider a binary system with two stars in a circular orbit. 
The stars are approximated as spherical masses Mi and M2 . 
Let a be the binary separation and 



G(Mi + M 2 ) 



n 1/2 



-1/2 



(2) 
(3) 



(1) 



the binary orbital frequency. 

Consider the binary frame rotating about the centre of 
mass with angular velocity ui, and let (R, <f>, z) be cylindrical 
polar coordinates such that star 1 (about which the disc 
orbits) is at the origin, while star 2 has fixed coordinates 
(a, 0,0). The axis of rotation passes through the centre of 
mass and is parallel to the z-axis. The Roche potential 

$ = -gm^r 2 + z 2 y 1/2 

-GM 2 (R 2 + a 2 - 2aRcos4> + z 2 ) 
- \u 2 R 2 + GM 2 a~ 2 R cos 

includes the gravitational potentials of the two stars and the 
centrifugal potential. 

The analysis of a tidally distorted disc is simplified 
by introducing non-orthogonal orbital coordinates (A, <f>, z) 
based on Paczyhski's orbits, instead of cylindrical po- 
lar coordinates. The region of the binary plane filled by 
Paczyhski's orbits is covered by a coordinate system (\,(f)), 
where A is a quasi-radial coordinate that labels the orbits, 
and 4> is the azimuthal angle. This naturally accounts for the 
principal tidal distortion of the disc. In this paper, only the 
case of a prograde circumstellar disc is considered, although 
the generalization to retrograde or circumbinary discs is 
straightforward. 

There are many possible ways to label the orbits. In the 
case of eccentric discs (Ogilvie 2001) it was natural to use 
the semi-latus rectum, which measures the angular momen- 
tum of the orbits, and reduces to R in the limit of circular 
orbits. However, the angular momentum is not constant on 
Paczyhski's orbits. A simple alternative is to define A to be 
the value of R at (f> — 0. In the absence of star 2, A reduces 
to R. In general, there exists a function R(X,4>), satisfying 
R(X, 0) = A, that defines the shape of the orbits and must 
be determined numerically. 

To carry out the vector calculus required in the fluid- 
dynamical equations, one must obtain the metric coefficients 
and connection components (cf. Ogilvie 2001). The two- 
dimensional line element is 

ds 2 = dR 2 + R 2 d<f> 2 = (R x dX + R^ d<f>) 2 + R 2 d(f> 2 , (4) 

where the subscripts on 7? stand for partial derivatives of 
the function R(X,(j>). Thus the metric coefficients are 

gxx = R\, g\c/> = RxR<i>, g^ = R 2 + R%- (5) 

As expected, the square root of the metric determinant is 
equal to the Jacobian of the coordinate system: 



1/2 r 
9 = J 



9{x,y) _ d(x,y) d(R, 



d(X,4>) d(R,<j>)d(X, 
The inverse metric coefficients are 

R<p 



g xx = 



R 2 + Rl 



R 2 R\ ' 



9 X * = - 



R 2 Rx' 



RR\. 



9 = 



R 2 ' 



(6) 



(7) 



The components of the metric connection, given by 

P? c = \g ad {d b g cd + d c g db - d d g bc ), (8) 

are found to be 



Rxx 



R\<t> R<t> 



AA R x ' ^ Rx R ' 
(R 2 + 2i?| — RR0<f,) 



rfx = 0, 



RRx 
r 4> _ Rx 



'1R ± 
R 



(9) 
(10) 
(11) 



The coordinate system is trivially extended to three di- 
mensions by incorporating the vertical coordinate z. With 
the exception of g zz — g zz = 1, all metric coefficients 
and connection components involving z vanish. The three- 
dimensional Levi-Civita tensor (alternating tensor) is 



abc 7—1 



(12) 



where e abc denotes the usual permutation symbol with val- 
ues 0, ±1. 



3 MODELLING THE TURBULENT STRESS 

It is now widely accepted that the turbulent stress in most, 
if not all, accretion discs is magnetohydrodynamic (MHD) 
in origin, resulting from the non-linear development of the 
magnetorotational instability (Balbus & Hawley 1998). In 
MHD turbulence any fluid variable X exhibits fluctuations 
X' about a mean value (X), so that X = (X) + X' with 
(X 1 ) = 0. Correlations between fluctuating quantities result 
in a non-zero mean turbulent Maxwell stress tensor, 



- — I BjBj — —BuBudi 
47r\ 1 J 2 k k ■ 



and Reynolds stress tensor, 



it 1 / / / . / . / / \ 

pUiUj + p [UiUj + UiUj + UiUj) 



(13) 



(14) 



(In this section only, Cartesian tensor notation is used, to- 
gether with the summation convention.) In numerical simu- 
lations the magnetic contribution is found to dominate and 
therefore any model of the stress should aim primarily to 
reflect the dynamics of the Maxwell tensor. 

In a recent paper (Ogilvie 2001) I proposed a simple 
viscoelastic model of the turbulent stress in an accretion 
disc. This was put forward as a straightforward generaliza- 
tion of the standard viscous model, used extensively in ac- 
cretion disc theory, to incorporate the effect of a non-zero 
relaxation time for the turbulence. It was also motivated 
by certain formal and physical analogies between MHD and 
viscoelasticity. Here I present a slightly refined version that 
accounts more completely for the transfer of energy between 
the mean motion and the turbulence, to heat and then to 
radiation. 

The turbulent stress is written as Tij — \T kk 8ij, where 
Tij is a symmetric second-rank tensor that can be associated 
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approximately with (_Bj'_Bj)/47rJ^] The equation of motion is 
then 



p(d t Ui + UjdjUi) 



-pdiQ-diip+lT^+diTv, 



(15) 



where p is the density, m the velocity, $ the external gravi- 
tational potential and p the pressure. The tensor Tij is taken 
to satisfy the equation 

+ rVTij = p(diUj + djUi) — (2p h + \p)(d k u k )8ij, (16) 

where r is the relaxation time, p the effective (dynamic) 
viscosity and pt, the effective bulk viscosity. The operator 
T> is a convective derivative acting on second-rank tensors, 
and is defined by 



VT 



ij — (dt + Ukdk)Tij - TikdkUj — TjkdkUi + 2Tijd k Uk-(17) 



With the aid of the equation of mass conservation, 

d t p + di( P Ui)=0, (18) 

an equation for the kinetic and gravitational energy of the 
fluid is then obtained in the form 



d t [p{\u 2 + $)] +8i [p(\u 2 + $)i 
= -Uidi{p + ^Tjj) + u^jTij 



(19) 



provided that <!> is independent of t. The turbulent energy 
density is identified as hTu (which can be associated approx- 
imately with {B'iB'i) /8ir), and this satisfies the equation 

dt(hTu) + di(TjjUi — TijUj) = Uidi(\Tjj) - mdjTij 



--{\Tu + 3p b diUi), 



(20) 



obtained by taking the trace of equation (|l6|). The first two 
terms on the right-hand side appear with the opposite sign in 
the preceding equation, and therefore represent the transfer 
of energy between the mean motion and the turbulence. The 
last term represents the loss of turbulent energy to heat. 
Accordingly the thermal energy equation is written as 



pT(d t s + Uidis) = -{\Tu + 3p h diUi) — diFt, 



(21) 



where T is the temperature, s the specific entropy and Ft 
the radiative energy flux. The total energy then satisfies an 
equation of conservative form, 

8t [p(\u + $ + e) + ±Tu] 

+di [p(\u + $ + h)ui + Tjjui - TijUj + Fi\ = 0, (22) 

where e is the specific internal energy, which satisfies the 
thermodynamic relation de = Tds — pd(l/p), and h — e + 
p/p is the specific enthalpy. 

Consider the following two limits of the above system. 
In the limit r — > the stress is related instantaneously to 
the velocity gradient by 

— \Tkk8ij = p{diUj +djm) + (p h - ^p)(dkUk)8ij. (23) 

The equation of motion then reduces exactly to the com- 
pressible Navier-Stokes equation. In the limit r — > oo, how- 
ever, the tensor Tij satisfies the equation 

3 The association is not perfect because equation ( |l6| ) does not 
necessarily guarantee that Tij is a positive semi-definite tensor, 
as (B<B'.)/4tt should be. 



TJTij = 0, (24) 

which implies that Tij is 'frozen in' to the fluid. This is pre- 
cisely the equation satisfied by the tensor BiBj/4-K when the 
magnetic field Bt is 'frozen in' to the fluid according to the 
induction equation of ideal MHD. The system of equations 
therefore reduces exactly to ideal MHD in this limit. 

For intermediate values of r, the viscoelastic model re- 
produces both the viscous (dissipative) and elastic (mag- 
netic) characteristics of MHD turbulence. It is intended to 
give a fair representation of the dynamical response of the 
turbulent stress to a slowly varying large-scale velocity field 
as found in a tidally distorted disc. The use of a tensor Tij 
in this model, rather than a vector or scalar field, is essential 
in order to allow off-diagonal stress components. Apart from 
the radiative energy flux, which in this paper is treated in 
the diffusion approximation, the viscoelastic model for r > 
results in a hyperbolic system of equations and is therefore 
'causal'. It also satisfies the continuum mechanical principle 
of material frame indifference. 

The viscoelastic model introduces a new dimensionless 
parameter for the disc. The Weissenberg number We can be 
defined as the product of the relaxation time and the angu- 
lar velocity, and is expected to be of order unity, although 
it has never been measured directly. In the case of a circu- 
lar, Keplerian disc the viscoelastic model predicts the usual 



'viscous' stress component Tn<f, 



f/if2 plus an additional 



component T^ — ^ We pQ. If the effective viscosity is given 
by an alpha prescription, p = otp/Q., the quantity T^^ can 
be associated approximately with a toroidal magnetic field 
with magnetic pressure | We a times the gas pressure. The 
turbulent energy is stored in this component. The shear en- 
ergy is tapped by the stress component Tr$, is converted 
into turbulent energy stored in T^, and then goes into heat 
and radiation. 



4 BASIC EQUATIONS 

The basic equations governing a fluid disc in three dimen- 
sions are now expressed in the orbital coordinate system. 
The equation of mass conservation is 



(8t + u a V a )p = -pV a 1i a , 



(25) 



where now u a is the velocity relative to the rotating frame. 
The equation of motion is 



[(dt + u b Vb)u a + 2r] abc L0 b u c 



= -pV"* - X7 a (p + \T\) + V b T 



nO-b 



(26) 



where uo a is the angular velocity of the rotating frame (only 
lu z — uj being non-zero). 

As described above, the tensor T ab satisfies the equation 

T ab + t [(dt + u c V c )T ab - T ac V c u - T bc V c u a + 2T ab V c u c 
= p(V a u b + VV) - (2p h + lp)(V c u c )g ab , (27) 

which, like the induction equation in MHD, is not affected 
by the rotation of the frame. The energy equation for an 
ideal gas is 



1 



7-1 



(dt + u a Va)P 



7 



7-1 
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+ -(lT a a + 3 l i b V a u a 



(28) 



For the fluid variables, introduce the expansions 



where 7 is the adiabatic exponent. The radiative energy 
flux, in the Rosseland approximation for an optically thick 
medium, is given by 



p a = _ 



16aT" 
3«p 



V a T, 



(29) 



where a is the Stefan-Boltzmann constant and k the opacity. 
The equation of state of an ideal gas, 

kpT 



P : 



p m m H 



(30) 



is adopted, where k is Boltzmann's constant, /i m the mean 
molecular weight and ijih the mass of the hydrogen atom. 
The opacity is assumed to be of the Kramers form 

k = C K pT- 7/2 , (31) 

where C K is a constant. The generalization to other power- 
law opacity functions is straightforward (Ogilvie 2000, 
2001). 

The effective viscosity coefficients are assumed to be 
given by an alpha parametrization. The precise form of the 
alpha prescription relevant to a tidally distorted disc is to 
some extent debatable. It will be convenient to adopt the 
form 



p, = ap 



-1/2 



Mb = "bP 



( GMx 
\ A 3 



-1/2 



(32) 



where a and a b are the dimensionless shear and bulk viscos- 
ity parameters. In the limit of a circular disc in the absence of 
star 2, this prescription reduces to the usual one, p = ap/Cl, 
etc., where Cl is the orbital angular velocity in the inertial 
frame. Finally, the Weissenberg number is defined by 



We 



= r( 



GMi 
V A 3 



1/2 



(33) 



where r is the relaxation time. Again, this reduces to We = 
Qr in the limit of a circular disc. For simplicity, the four 
dimensionless parameters of the disc, (a, ah, We, 7), will be 
assumed to be constant throughout this paper. The other 
parameter in the problem is the mass ratio q = M2 /Mi . 



5 SOLUTION OF THE EQUATIONS 

The analysis proceeds on the basis that the disc is almost 
steady in the binary frame, and evolves only on a time-scale 
(the viscous time-scale) that is long compared to the orbital 
time-scale. It is also assumed that the fluid variables vary 
smoothly in the horizontal directions, not on a length-scale 
as short as the semi-thickness of the disc. Under these con- 
ditions, the temporal and spatial separation of scales allows 
the fluid dynamical equations to be solved by asymptotic 
methods. 

Let the small parameter e be a characteristic value of 
the angular semi-thickness H/R of the disc. Adopt a system 
of units such that the radius of the disc and the character- 
istic orbital time-scale are 0(1). Then define the stretched 
vertical coordinate ( = z/e, which is O(l) inside the disc. 
The slow evolution is captured by a slow time coordinate 
T = e 2 t. 



u = 


e 2 u^(X,<t>,(,T) + 0(e 4 ), 


(34) 


„> - 

u — 


ll(A, <p) + e u 2 (A, <p,(,,l ) + U(e ), 


(35) 


z 

u = 


euU\,d>,C,T) + Ofe 3 ), 


(36) 


p = 


e s [po(A,<?i,C,T) + 0(e 2 )] , 


(37) 


V = 


e s+2 [p (A>,C,T) + O(6 2 )] , 


(38) 




e s+2 J To ^( Aj ^ CjT ) + ( e 2)] f 


(39) 


rj~)\(p 


e s+2 [ ro ^(A>,(,T) + 0(e 2 )] , 


(40) 


rp4>4> _ 


s+2 \ T t<t>tx,<t>,C,T) + 0(e 2 )\ , 


(41) 




e s+3 [T 1 Az (A,0,C,T) + O( £ 2 )] , 


(42) 




e s+3 [Tt(\^,C,T) + 0(e 2 )] , 


(43) 




e s+2 [r "(A,^,C,T) + O( £ 2 )] , 


(44) 


H = 


e s+2 [ MO (A,0,C,T) + O(6 2 )] , 


(45) 


Mb 


e s+2 [ Mbo (A,0,C,T) + O( e 2 )] , 


(46) 


7 = 


e 2 [To(A,^,C,T)+0(6 2 )] , 


(47) 


F z = 


e s+3 [F (X,^C,T)+O(e 2 )] , 


(48) 



while the horizontal components of F are 0(e s+4 ). Here s 
is a positive parameter, which drops out of the analysis, 
although formally one requires s = 11/3 in order to bal- 
ance powers of e in the opacity law. Note that the dominant 
motion is an orbital motion with angular velocity £1 in the 
rotating frame, independent of £ and T. 

The Roche potential is expanded in a Taylor series 
about the mid-plane, 



$ = $o(A,<W + i e 2 C 2 $ 2 (A,0) + C>(e 4 ), 
where 

$0 = -GM 1 R~ 1 - GM 2 (R 2 + a 2 - 2aRcos(j,y 1/2 



-ui R + GM2C1 Rcosc 



$ 2 = GM1R- 3 + GM 2 {R 2 + a 2 - 2a.Rcosc 



-3/2 



(49) 

(50) 
(51) 



When these expansions are substituted into the equa- 
tions of Section ^, various equations are obtained at different 
orders in e. The required equations comprise three sets, rep- 
resenting three different physical problems, and these will 
be considered in turn. 



5.1 Orbital motion 

The horizontal components of the equation of motion (|26| 
at leading order [0(e s )] are 



Po (r^n 2 - -^cun^j = -p (g xx d^o + g^d^o), (52) 

p (nd*n + r^n 2 + ^wn) = -po(s a ^a<i>o+<7^$o).(53) 

As $0 is naturally given as a function of R and <j>, these may 
be written in the form 



(R 2 + 2R% - RR^)tt 2 + 2(R 2 + J?;)wfi : 

,(d$ \ R4, f dV 
V dR J 4, R \ d</> , 



(54) 
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Ud^ [R 2 (n + uj)] 



(55) 



Since no derivatives with respect to A appear, these equa- 
tions may be solved as a third-order system of ordinary dif- 
ferential equations (ODEs) on < 4> < 2n for each orbit 
separately. It is convenient to non-dimensionalize the equa- 
tions by taking a as the unit of length and (GAfi/o 3 ) 1 ' 2 as 
the unit of angular velocity. The unknown functions R and 
il are subject to periodic boundary conditions. The Jacobi 
constant, 



£=i(ii 2 +i4)fi 2 + $o, 



(56) 



is an integral of the system, satisfying E$ = 0. 

The solution of these equations amounts to solving for 
Paczyhski's orbits, and yields quantities such as R, R$, R^, 
Q and c^fi. By differentiating the equations with respect 
to A, one obtains a further set of equations that determine 
quantities such as R\, R\^ and d\fl, which will all be re- 
quired below. 



5.1.1 Surface density 

The equation of mass conservation (^^) at leading order 
[0(e")] is 



(0.84, + uld^po = -p 



J 



d<f,(JFl) + d^u\ 



(57) 



Introduce the surface density at leading order [0(e s+1 )], 
E(A,</>,T) = / podC (58) 



Then the vertically integrated version of equation (^) may 
be written in the form 



d^jtQ) = 0, 



(59) 



which determines the variation of surface density around 
the orbit due to the tidal distortion. It will be convenient to 
introduce a pseudo-circular surface density, 



S(A,T) = — / J£d</>, 



(60) 



which has the property that the mass contained between 
orbits Ai and A2 is 

2tt / SAdA. (61) 

J Ai 

5.1.2 Kinematic quantities 

The shear tensor of the orbital motion will be required be- 
low. This is defined in general by 



S 



ab 1 /yjo, b 



( W + v" 



and has the expansion 
S ab = Sg h +0(e). 



(62) 
(63) 



4 Throughout, integrals with respect to f are carried out over 
the full vertical extent of the disc, and integrals with respect to 
4> from to 2iz. 



The horizontal components at leading order are 

Q xx {R 3 R\4> + RR\(t>R% + R\R% — RRxR^R^O. 
60 = Wr\ ' 



(64) 



~,\$ (-R 2 + R\ 
5 ° ~ 2R?R\ 



Re 



+ 



(RRxRcf}^ — RR<j>R\4> 



2R 2 R X 

2R\R^)Q. 



°0 — 



2R 3 R\ 



The orbital contribution to the divergence \7 a u a is 



A = -<9 (Jfi) 



~ \ A3 ) 



gi, 



(65) 
(66) 

(67) 



where g\ (<j>; A) is a dimensionless quantity, equal to zero for 
a circular disc, and which can be evaluated numerically from 



the solution of the problem in Section 5.1 



Note, however, that there is also a vertical shear com- 
ponent Sq z — d(uf of the same order, which remains to be 
evaluated. This contributes to both the divergence and the 
dissipation rate at leading order. 

Let 



/GMA 1 ' 2 
Q = I — r=— I gn, 



A 3 

GMi n 1/2 



(68) 



(69) 



where gn(<f>;X) and <7u(A) are dimensionless quantities. In 
the limit of a circular disc, gn + g w = 1. 
Noting that 



dej) 

is the orbital period, one may write 



(70) 



(71) 



where 32(0; A) is a dimensionless quantity, equal to unity for 
a circular disc, and satisfying 



2ttA 
92 = JOP' 

gadding? = -51. 

Further useful quantities are defined according to 
A^O=f^V /2 



V A 3 ' 53 ' 



1/2 



/ GMi\ 

$ 2 = (^r) g, 



9-t- 



and 



(GMi\ 



1/2 



5T = {^) 
so that (gs, g<t, g$, s AA , As A< *, A 2 s*^) are dimensionless. 



(72) 
(73) 

(74) 
(75) 
(76) 

(77) 
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5.2 Vertical structure and vertical motion 

The equation of mass conservation (^H|) at leading order 
[G(e s )] is 



(fic^ + u{d c )p = -po(A + d ( uf). 
The energy equation (E8I) at leading order [0(e 



(78) 



3 + 2\ 



7 — 11 " I 7- 1 



po(A + d ( ut) 



+6fj, b0 (A + d ( ui)]-d ( F . 



(79) 



The vertical component of the equation of motion ( |26[ ) at 
leading order [0(e s+1 )] is 

po(Q.d<t, + uld c )u\ = -p <I>2C 

-8 C [p + i (gxxT x x + 2g^T x * + g u T$* - I?*)](80) 

The required components of the stress equation (E?]) at lead- 
ing order [0(e s+2 )] are 

t xx + r [(nfy + «;a c )2^* + 2(A + Scu*)^] 

= 2 M S„ AA - (2/i b0 + |/io)(A + d c ul)g xx , (81) 

t a,?> + r [(nfy + u;e f )2£* - r, AA a A n - 

+2(A + 9 C ^)T A0 ] 
= 2im>SS* - (2/i b0 + |/io)(A + c «")ff A *. (82) 

2^* + r [(flfy + uld<)T*+ - 2T A ^ A fi - zz*%n 
+2(A + C «?)2?*1 



= 2f*>Sj* - (2Mb.. + f Po)(A + 9 <Ul ) 5 

T<f + r [(fify + u?0 f )T?" + 2AT "] 
= 2p d(ul - (2/ibo + §Mo)(A + 
The constitutive relations at leading order are 



F = - 



16aT f 



13/2 



(1 



PO 



kp T 



-d c T , 



(83) 
(84) 

(85) 
(86) 



po = «Po 



p m m H 

fGMi\~ 1/2 /GMi\~V a 
^^3- J - W>o = «bPo (^3- J • (87) 

The solution of equations (frs])-(^7^ is found by a non- 
linear separation of variables, a similar method to that used 
for warped discs (Ogilvie 2000) and eccentric discs (Ogilvie 
2001). As an intermediate step, one proposes that the solu- 
tion should satisfy the generalized vertical equilibrium rela- 
tions 

d Po , / GAfi \ 



dF . 9a /GMi\ 

~9c = fl T Po [-J^) 



1/2 



(89) 



in addition to equations (|85[)-(|87j). Here fi(<f>; A) and /2(</>; A) 
are dimensionless functions to be determined, and which are 
equal to unity for a circular disc with We — * 0. Physically, 
fx differs from unity in a tidally distorted disc because of 



the enhanced dissipation of energy and because of compres- 
sive heating and cooling (equation |79|) . Similarly, /2 reflects 
changes in the usual hydrostatic vertical equilibrium result- 
ing from the vertical velocity, the turbulent stress and the 
variation of the vertical oscillation frequency around the or- 
bit (equation ^) . Following Ogilvie (2000) , one identifies a 
natural physical unit for the thickness of the disc, 



^9o^ x / 14 E 3 /14 ^GMij " 3/7 - 



15/28 



and for the other variables according to 



U P = St/" 1 , 



/ GM\ \ rr2 



V A3 



Uf = 



/9a \ fGMi \ 1/2 
V 4 ) { A3 



UhU p 



(90) 

(91) 
(92) 

(93) 
(94) 



The solution of the generalized vertical equilibrium equa- 
tions is then 



r-t fl/14 f -13/28 3/14,-, 
C — £* Jl h 92 U H, 

If \ r-1/14 *13/28 ll/14 rr 

Po = P»(C*)/i Si 9 2 U p , 

rf . \ ,1/14 „15/28 17/14 Tr 
P0 =P»(C*)A S 2 92 U p , 

rr\ rr I j- \ rl/7 fl/14 3/7 rr 

To = r„(c*)/i' f 2 9 2 u T , 

Ei t/- ^ /8/7.1/14 10/7 r » 

To = F,(C*)/i / 2 5 2 ^f, 



(95) 
(96) 
(97) 
(98) 
(99) 



where the starred variables satisfy the dimensionless ODEs 

(100) 

(101) 



dp, 

dc: = - 



dF» 



dC* 

p* = p,T,, 
subject to the boundary conditions 
F,(0) =p.(C B .) =T.(Cs.) =0 
and the normalization of surface density, 

PCs. 



p. dC* = 1. 



(102) 
(103) 

(104) 
(105) 



Cs. 



Here is the dimensionless height of the upper surface of 
the disc. The solution of these dimensionless ODEs is easily 
obtained numerically in a once-for-all calculation (Ogilvie 
2000), and one finds £ s * ~ 2.543. However, the details of the 
solution do not affect the present analysis.^] 

One further proposes that the vertical velocity is of the 

form 



J It is assumed here that the disc is highly optically thick so 
that 'zero boundary conditions' are adequate. The corrections for 
finite optical depth are described by Ogilvie (2000). 
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. . / GMi \ 1//2 

%=/3 ( — J 



c, 



(106) 



where /3(<^>;A) is a third dimensionless function, equal to 
zero for a circular disc, and that the stress components at 
leading order are of the form 

TS b =t ab (4>;X)p , (107) 

so that (t xx ,\t X4 ',\ 2 t 4 " l ',t zz ) are dimensionless coefficients. 

When these tentative solutions are substituted into 
equations (|7^)-(|84|) one obtains a number of dimensionless 
ODEs in 4>, which must be satisfied if the solution is to be 
valid. From equation (JT^) one obtains 

gn (-d^hxfx + yd ln/ 2 ) = -14/ 3 - 3g t . (108) 

Similarly, from equation ( pFs] ) one obtains 

+T nXT [d xxtXX + 2 9^ tM + + t zz + 6ab(ffi + h)] 



2We 
9 f 



(109) 



From equation (g0|) one obtains 

gnd^h = -fi ~ 35 

+h [1 + \{gxxt xx + 2g M t x * + 9H ^ - t zz )] . (110) 

Finally, from equations (jsi|)-([s4|) one obtains 

t xx + We [<7n(c^i AA + t xx d4, In f 2 ) + (3 f 3 + gi)t xx ] 

= 2as xx - (2a b + §a)Cffi + / 3 ) 3 AA , (111) 

\t X4> + We [<?n(^(At A0 ) + A* A % In /a) 

+(3/a + gi )\t x + - g s t xx - 9i \t x *] 
= 2aXs X4, - (2a h + |a)(ffi + / 3 )A ff A ^, (112) 

A a t** + We [fln(S*(A a t**) + AV%ln/ 2 ) 

+(3/a +gi)\ 2 t^ - 2g 3 \t x * - 2g,X 2 t^] 
= 2aA 2 s^ - (2a b + |a)( fl i + / 3 )A a g* (113) 



t zz + We [ffn(^t zz + t zz d^ In / 2 ) + (/ 3 + gi )t zz ] 
= 2a/ 3 - (2ab + |a)(gi + /s). 



(114) 



These ODEs should be solved numerically for the 
functions (/i, / 2 , / 3 , t AA , A^, A 2 ^ , t zz ) subject to periodic 
boundary conditions /i(27r;A) = /i(0;A), etc. Note that, 
in the limit of a circular disc (gi = 0, <? 2 = 1, g 3 = — §, 



<?4_= 0, <?5 = 1, g AA _ = 1, A ff A ^ = 0, AV* = 1, s\\ = 0, 
A sxtf, = — |, A 2 s^0 = 0), the solution is /i = 1, 
/ 2 = (1 + fWea)" 1 , / a = 0, t AA = 0, A^ = 



-fa, 



5.3 Slow velocities and time-evolution 

It is possible to expand further in order to determine the 
mean accretion flow in the disc and deduce the equation 
governing the evolution of the pseudo-circular surface den- 
sity E(A,T). This is a generalization of the diffusion equa- 
tion for the surface density of a standard, circular accretion 



disc, which includes the effects of tidal distortion and tidal 
torques. The tidal truncation of the disc should be discussed 
in the context of this evolutionary equation. 

The method is again similar to that used for eccentric 
discs (Ogilvie 2001), but is complicated by the fact that the 
specific angular momentum is not constant on Paczyhski's 
orbits. As this aspect of the problem introduces much addi- 
tional complexity without affecting the main conclusions of 
this paper, it is deferred for future work. 



6 INVESTIGATION OF THE SOLUTIONS 
6.1 General remarks 

The preceding analysis is implemented numerically as fol- 
lows. First, one selects the binary mass ratio q and the di- 
mensionless disc parameters (a, otb, We, 7).^] A discrete set 
of Paczyhski's orbits is computed, with equally spaced val- 
ues of A from the inner radius of the disc to the firs t poin t 
of orbital intersection. The dimensionless ODEs (108)— (114) 
are then solved to yield the functions fi for each orbit. 

The solutions may be interpreted by considering phys- 
ical quantities of particular interest, such as the surface 
brightness of the disc and its semi-thickness. The surface 
brightness (in continuum radiation) is the value of the ver- 
tical radiative flux F z at the upper surface. From equation 
( p9[ ) this may be written in the form 

F+ = f F C F a^ 7 (£|*L) 9/14 E 10//7 , (115) 
where 

O.- (»)«%.«-> <ne> 

is a constant depending only on the opacity law, and 
/f(0; A) is a dimensionless function defined by 



f _ f8/7, 1/14 10/7 

If — Ji J 2 3 2 1 



(117) 



and equal to unity for a circular disc with We — > 0. This con- 
tains the azimuthal variation of the surface brightness and 
represents the principal tidal distortion of the disc. It does 
not depend on £(A, T). The additional variation of F + with 
A an d T is contained in the remaining factors in equation 
(115), and requires a knowledge of E(A,T). 
Similarly, the semi-thickness of the disc is 



where 



C H 



(~3C, 



(118) 



(119) 



is a constant depending only on the opacity law, and 
fn {4>\ A) i s a dimensionless function defined by 



f _ /1/14 .-13/38 3/14 
JH — ,/l h 92 > 



(120) 



6 There is no difficulty in principle in extending the analysis of 
this paper to situations in which these dimensionless parameters 
vary within the disc. 
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and equal to unity for a circular disc with We — > 0. 

There are therefore two aspects to a complete solution 
for a tidally distorted disc. The first aspect concerns the de- 
termination of Paczyriski's orbits and the functions fi, which 
indicate how the physical quantities vary around each orbit. 
This part of the problem does not depend on the distri- 
bution £(A, T) of surface density over the set of orbits. In 
this sense, the non-axisymmetric tidal distortion of the disc 
is fixed and depends only on g and the dimensionless disc 
parameters (a, at,, We, 7). However, to obtain the correct 
variation of physical quantities from one orbit to the next, 
and in time, requires a knowledge of E(A,T). This second 
aspect of the problem requires a solution of the evolutionary 
equation described briefly in Section 5.S . During the course 
of a dwarf- nova outburst, the first aspect of the problem 
remains fixed except inasmuch as a varies, but the second 
aspect changes significantly as the surface density evolves. 

6.2 Non-axisymmetric structure 

The non-axisymmetric structure of a tidally distorted disc 
may be illustrated by plotting the functions /f(</>;A) and 
/#(</!>; A). This is done in Fig. 2 for a binary mass ratio q = 
0.5 and illustrative disc parameters a — 0.1, ax, = 0, We = 
0.5, and 7 = 5/3. This represents a disc with a Kramers 
opacity law and a moderate viscosity and relaxation time. 

These results are very typical. Both the surface bright- 
ness and the thickness display a dominant m — 2 distor- 
tion in the outer part of the disc. The phases are different, 
however. While the brightened sectors of the outer disc are 
approximately aligned with the companion, the thickened 
sectors are misaligned by approximately 45°. When trans- 
lated into velocity space, the brightened sectors have the 
wrong phase to explain the tomograms of IP Peg and other 
systems, while the thickened parts have the correct phase. 

6.3 Emission lines and synthetic Doppler 
tomograms 

This suggests that the observed patterns can be explained 
if the line emission is correlated with the thickness of the 
disc rather than the continuum brightness. This is consistent 
with the notion that the emission lines originate in parts of 
the disc that are elevated and therefore irradiated from the 
centre of the system by the white dwarf, boundary layer 
and/or inner disc (Robinson et al. 1993). 

Simple synthetic Doppler tomograms based on this 
premise can then be constructed by plotting an image rep- 
resenting the excess elevation of the disc in velocity space. 
The quantity plotted in Fig. 3 is fn — }ho, where 

/ O \ 13/28 

/ho= (l+-Weaj 

is the constant value of fjj in the absence of star 2. The 
velocity field used in constructing the tomograms is the or- 
bital motion of Paczyriski's orbits, translated into the iner- 
tial frame. 

These synthetic tomograms should not be considered as 
a serious attempt to reconstruct the observed data, such as 
that made by Steeghs & Stehle (1999). This would require a 
better understanding of the mechanism of line emission and 
a proper treatment of the instrumental and data-analytical 



limitations of Doppler tomography. Rather, Fig. 3 should be 
taken as an indication of the probable locations in velocity 
space of line emission under the irradiation hypothesis. 

The dependence of the two-armed patterns on the disc 
parameters is rather subtle. Increasing the values of a, o?b or 
We can cause phase shifts that rotate the patterns to some 
extent. Decreasing 7 moves the vertical resonance outwards 
in real space (inwards in velocity space) and can also increase 
the amplitude of the distortion by making the disc more 
compressible. 



6.4 The role of the vertical resonance 

The significance of the vertical resonance in this problem can 
be demonstrated by considering an approximate model that 
corresponds to a linear, inviscid analysis of the full equa- 
tions. In the limit a, — > the dimensionless turbulent 
stress coefficients t ab tend to zero and one obtains the much 
simpler system of ODEs 



gn (-84, In /1 + y 84, In / 2 ) = -14/ 3 - 3fli, 
'7+1" 



gn 



d<t> m h — 



/s - 91, 



-fl 



+ h 



(122) 

(123) 
(124) 



When the tidal distortion is weak, one may write / 2 = 1 + 
0/2, gn = (l~g^)+6gn and g 5 = l + 5g 5 , where 5f 2 , f 3 , Sgn, 
gi and Sgs are small quantities. In the linear approximation, 
products of small quantities may be neglected, and this leads 
to a linear, inhomogeneous equation of second order for 8/2, 



[(l-.9c) 2 a| + ( 7 + l)] Sf 2 =X, 
where 

X = (7 + l)5 g5 - (7 - 1)(1 - g u )d^ gi . 



(125) 



(126) 



The forcing function may be resolved into Fourier compo- 
nents, 



X(0;A) = Re^X m (A)e im0 , 



(127) 



where the coefficients X m are real by virtue of the symmetry 
of Paczyriski's orbits. The solution is then 



00 



5/ 2 = Re ^ [( 7 + i)- m 2 (i-^) 2 ] L X m e im t (128) 

m=0 

It is singular at the vertical resonances, where the quantity 
in square brackets vanishes. Now 



(121) ^ = (i + g) 1/2 (^) 



M 3/2 



(129) 



and so the inner/outer resonance for mode m > occurs at 



H 1 

a 



■q) 



-1/3 



IT -(7 + 1) 



1/2 



2/3 



(130) 



as determined by Lubow (1981). In particular, for q = 0.5 
and 7 = 5/3, the m = 2 inner vertical resonance occurs at 
A/a » 0.282. 

The phase of the response to m = 2 forcing changes by 
n as one passes through the radius of the resonance. In the 
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Figure 2. Non-axisymmetric distortion of a disc with parameters a = 0.1, «b = 0, We = 0.5 and 7 = 5/3 in a binary with q = 0.5. 
The axes are in units of the binary separation, with the companion located at (1, 0). The thick line represents the outer edge of the disc. 
Left: contours of the flux distortion fp, with spacing 0.2. The dotted, dashed and solid contours represent values less than 1, equal to 1 
and greater than 1, respectively. Right: contours of the height distortion /jy, with spacing 0.05. 



linear, inviscid limit, the response is formally singular at the 
resonance. In reality, the response is limited by the emission 
of a wave (Lubow 1981), or by non-linearity or dissipation 
(Ogilvie, in preparation). The singularity is then avoided 
and a smooth variation of amplitude and phase occurs. 

The distortion of the thickness of the disc is sensitive 



mainly to the quantity /2, according to equation (120) 



Therefore the vertical distortion is amplified in the vicin- 
ity of the vertical resonance and undergoes a phase shift as 
one passes through the resonance. This effect can be seen 
in Fig. 2 (right panel). For the m = 2 component of the 
solution, the phase shift of tv corresponds to a rotation of 
the peak through 90°. This lends a certain 'spirality' to the 
pattern in the vicinity of the resonance. However, the solu- 
tion does not represent a wave or shock, and the shape of 
the pattern is independent of the thickness of the disc. 



7 SUMMARY AND DISCUSSION 

In this paper I have constructed a detailed semi-analytical 
model of the non-linear tidal distortion of a thin, three- 
dimensional accretion disc by a binary companion on a cir- 
cular orbit. The analysis allows for vertical motion and ra- 
diative energy transport, and introduces a simple model for 
the turbulent magnetic stress. The fluid dynamical equations 
are solved by means of a consistent asymptotic expansion, 
using the fact that the disc is thin. The solution is formally 
exact in the limit of a thin disc. 

The tidal distortion affects quantities such as the sur- 
face brightness of the disc and its thickness. As expected, 
the distortion is greatest in the outer part of the disc and 
is predominantly two-armed. An important feature of the 
solution is the m = 2 inner vertical resonance, first anal- 
ysed by Lubow (1981), which typically lies within the tidal 



truncation radius in binary stars. This intrinsically three- 
dimensional effect has a significant influence on the ampli- 
tude and phase of the tidal distortion. 

I have shown that the two-armed features observed in 
Doppler tomograms of the discs of dwarf novae in outburst 
find a natural explanation if the emission is associated with 
the tidally thickened sectors of the outer disc that are ele- 
vated and may therefore be irradiated by the white dwarf, 
boundary layer and/or inner disc during the outburst. There 
is a slight spirality to the structure, arising from the phase 
shift that occurs on passing through the vertical resonance. 
There is some dependence of the pattern on the mass ratio 
and the dimensionless parameters (a,cib,We, 7) of the disc 
(defined in Section [|). In particular, the adiabatic exponent 
7 determines the location of the vertical resonance. Future 
observations of higher quality might help to constrain the 
values of these parameters. In principle, the distortion of 
the continuum surface brightness of the disc, which displays 
a different pattern (e.g. Fig. 2), might also be seen in eclipse- 
mapping observations. 

This possible interpretation of the observed features is 
fundamentally different from one based on spiral waves or 
shocks. There are no waves or shocks in the present analysis. 
The pattern of emission merely reflects the tidal distortion 
of the disc and is independent of the angular semi-thickness 
H/r, in the limit that this is small. The observability of 
the pattern will nevertheless vary during the outburst cycle 
depending on the spreading of mass to the outer parts of the 
disc and, in the irradiation hypothesis, on the luminosity of 
the central regions. 

In a recent paper, Smak (2001) has presented argu- 
ments against the interpretation of the observed features 
as spiral shocks. He argues in favour of the irradiation hy- 
pothesis and suggests that sectors of the outer disc may be 
elevated through the horizontal convergence of the three- 
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Figure 3. Simple synthetic tomograms for a disc in a binary. The velocity scale is based on the orbital parameters of IP Peg. The plotted 
intensity is proportional to the elevation of the disc above its undistorted value. Top left: standard parameters, as in Fig. 2. Top right: 
standard except that a = 0.3. Middle left: standard except that We = 0.1. Middle right: standard except that We = 1. Bottom left: 
standard except that 7 = 1.2. Bottom right: standard except that q = 0.1 (the velocity scale assumes the same value of GM\/a as in IP 
Peg). 
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body orbits. The present paper supports the arguments of 
Smak (2001) to a certain extent, but within the context of 
a three-dimensional fluid dynamical model. The principal 
factors lead ing to the thickening of the disc can be seen in 
Section 6.4 . Vertical expansions and contractions are driven 



by the variation of the vertical gravitational acceleration 
around the tidally distorted orbits and, to some extent, by 
the horizontal convergence of the flow in the binary plane. 
The response of the disc to this forcing is regulated by the 
inner vertical resonance. 

The analysis in this paper could be extended and im- 
proved in a number of ways. Mass input to the disc, its 
ultimate tidal truncation, and the spreading of mass dur- 
ing an outburst, should properly be discussed within the 
context of the evoluti ona ry equation for the surface density, 
alluded to in Section 5.3 but not derived here. The solution 



of this equation would determine the relative weighting to 
be assigned to the different orbits. The viscoelastic model 
for the turbulent stress, although an improvement on purely 
viscous treatments, could also be refined. Finally, a better 
understanding of the mechanism of line emission would en- 
able a closer connection to be made between the computed 
models and the observational data. 
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